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The  detailed  description  of  the  vibrational  energy  content  is  important  for  many  gas 
flow  processes  including  experimental  gas  measurement  techniques,  shock  layer  vibration- 
dissociation  coupling,  and  vibrational  energy  freezing  in  strong  expansions.  Various  vibra¬ 
tional  relaxation  models  have  been  developed  for  use  within  the  direct  simulation  Monte 
Carlo  (DSMC)  method.  This  study  compares  the  Larsen-Borgnakke  and  Forced  Harmonic 
Oscillator  vibrational  relaxation  models  through  comparison  with  detailed  quasi-classical 
trajectory  (QCT)  results  for  an  atom-diatom  system  with  focus  on  the  effects  of  the  form 
of  the  relaxation  model.  These  QCT  results  have  been  extensively  compared  with  mea¬ 
sured  data  and  provide  a  unique  opportunity  to  evaluate  these  inelastic  models.  This  work 
also  compares  the  results  from  each  model  for  a  vibrationally  exciting  relaxation  problem. 
Even  when  the  two  models  have  consistent  vibrational  collision  numbers  to  match  bulk 
temperature  predictions  for  relaxation  problems,  the  predicted  evolution  of  the  vibrational 
distribution  function  differs.  This  disagreement  results  in  different  prediction  of  radiative 
emission  spectra. 


I.  Introduction 

Excitation  and  quenching  of  vibrational  energy  modes  through  collision  relaxation  is  an  important  pro¬ 
cess  in  a  wide  range  of  high  speed  flows.  High  temperature  shock  layers  that  form  in  front  of  hypersonic 
vehicles  can  lead  to  significant  vibrational  excitation.  Due  to  the  relatively  slow  vibrational  excitation  rates, 
the  vibrational  energy  may  only  be  partially  excited  compared  to  the  translational  and  rotational  energy 
modes  within  most  of  the  shock  layer  and  the  vibrational  energy  distribution  function  (VEDF)  may  not  be 
characterized  by  an  average  mode  temperature  due  to  the  strong  non-equilibrium  effects.  In  high  expansion 
flows,  the  vibrational  energy  may  remain  frozen  at  near-stagnation  conditions  and  large  regions  of  nozzle 
flows  may  experience  highly  non-equilibrium  VEDFs.  This  non-equilibrium  VEDF  has  a  strong  effect  on 
the  coupled  vibration-dissociation  of  molecular  species.  In  addition,  the  non-equilibrium  VEDF  has  a  strong 
effect  on  measured  radiative  emission  that  is  used  to  estimate  vibrational  energy  content  in  experimental 
systems.  The  non-equilibrium  VEDF  is  also  important  to  evaluate  fundamental  gas  dynamic  properties, 
such  as  the  ratio  of  specific  heats  and  thermal  conductivity. 

Many  vibrational  relaxation  models  have  been  developed  to  use  in  both  continuum  flows  for  use  in 
traditional  Computational  Fluid  Dynamics  (CFD)  and  non-continuum  flows  for  use  with  rarefied  flow  de¬ 
scriptions,  such  as  the  direct  simulation  Monte  Carlo  (DSMC)  method.  At  non-continuum  (or  rarefied) 
conditions,  the  non-equilibrium  effects  of  vibrational  relaxation  become  even  more  important  as  there  are 
insufficient  collisions  to  equilibrate  the  fast  relaxation  of  the  velocity  distribution  function  (VDF)  as  well  as 
the  relatively  slow  relaxation  of  the  VEDF.  Recently  Boyd  and  Josyula,1  have  compared  standard  DSMC 
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vibrational  relaxation  approaches2,3  with  a  recently  introduced  high-fidelity  vibrational  relaxation  model.4 
They  also  have  compared  the  high-fidelity  results  with  thermally-averaged  results  obtained  by  Billing  and 
Fisher.5 

In  this  work,  we  extend  the  previous  work1  by  comparing  both  low-  and  high-fidelity  results  to  state- 
resolved  quasi-classical  trajectory  (QCT)  results  from  Braunstein  and  Duff6.  These  results  have  been  ex¬ 
tensively  compared  with  experimental  measurements7  and  provide  a  unique  opportunity  to  evaluate  a  wide 
range  of  vibrational  inelastic  collision  models.  First,  we  present  the  ro-vibrationally  resolved  cross  sections 
of  Braunstein  and  Duff.  Then,  we  describe  the  low-  and  high-fidelity  models  that  are  used  within  the  DSMC 
method  with  focus  on  the  inherent  form  that  dictates  the  VEDF  that  is  sampled  in  vibrationally  inelastic  col¬ 
lisions.  Next,  we  compare  the  models  to  available  data  from  the  QCT  modeling  for  vibrational-translational 
energy  exchange  for  an  atom-diatom  system.  In  particular,  we  assess  the  models  by  comparing  the  predicted 
VEDF  during  a  vibrational  inelastic  collision  at  a  given  initial  state  and  collision  velocity.  Finally,  the  two 
models  are  applied  in  the  DSMC  method  to  predict  vibrationally  excitation  in  an  adiabatic  relaxation  case, 
where  both  models  are  tuned  to  agree  in  the  bulk  relaxation  process.  This  isolates  the  effect  of  the  model 
form  on  the  predicted  evolution  of  the  VEDF  and  associated  radiative  emission  spectra.  Detailed  comparison 
of  the  evolution  of  the  VEDF  and  emission  spectra  are  compared  for  each  vibrational  relaxation  model.  The 
paper  closes  with  summary  and  conclusions. 

II.  Available  State-to-State  Transition  Data 

In  this  section,  we  present  model  results  from  quasi-classical  trajectory  (QCT)  calculations  of  Braunstein 
and  Duff.6  These  results  have  been  extensively  compared  with  experimental  measurements7  and  provide 
an  opportunity  to  evaluate  a  wide  range  of  vibrational  relaxation  models.  We  present  the  ro-vibrational 
cross  sections  for  a  number  of  relative  collision  velocities  and  compare  the  post  vibrationally-inelastic  VEDF 
(PVI-VEDF),  which  is  described  in  Section  III,  at  selected  collision  classes. 

In  Braunstein  and  Duff,6  the  potential  energy  surfaces  of  the  three  lowest  electronic  triplet  states  of 
CO2  which  lead  to  0(3P)  +  CO(1E+),  3  A/,  13A",  and  2 3A",  were  computed  at  the  complete-  active-space- 
self-consistent-field  plus  second-order  perturbation  theory  (CASSCF-MP2)  level  with  a  modest  631  +  G(d) 
basis.  Potential  energy  surfaces  were  fit  with  a  global  functional  form.  These  potential  energy  surfaces 
were  used  to  run  classical  trajectory  calculations  for  ro-vibrationally  resolved  0(3P)  +  CO(1E+)  collisions. 
When  O  +  CO  collide,  there  is  direct  collisional  excitation  of  the  CO  and  there  is  also  reaction  where  the 
incoming  O  exchanges  with  the  O  in  the  reactant  CO.  Both  collisional  excitation  and  reaction  channels  lead 
to  vibrational  excitation  of  the  product  CO.  The  QCT  results  presented  below  include  both  reactive  and  non¬ 
reactive  channels  summed  together.  Their  results  agree  well  with  available  angle  and  energy  resolved  cross 
section  measurements7  except  at  low  temperature.  Comparisons  have  been  made  with  measured  vibrational 
excitation  cross  sections  and  infrared  emission  spectra  of  the  nascent  CO  products  at  3.4  eV  collision  energy. 

It  is  important  to  note  that  the  methods  used  to  generate  these  results  are  unable  to  accurately  predict 
the  vibrationally  elastic  collision  cross-section.  In  this  study,  a  separate  total  cross  section  model  is  used  in 
conjunction  with  the  QCT  results  to  form  an  overall  relaxation  rate.  In  this  study,  this  is  not  required  where 
the  focus  is  on  the  PVI-VEDF.  These  results  provide  an  opportunity  to  evaluate  a  wide  range  of  inelastic 
vibration  models. 

Figure  1  shows  the  post  vibrationally-inelastic  ro-vibrational  energy  distribution  function  (RVEDF)  while 
varying  the  initial  vibrational  level  and  rotational  temperature  of  the  CO  molecule  at  two  relative  velocities. 
Figures  1(a),  (c),  and  (e)  show  the  post  vibrationally-inelastic  RVEDF  at  moderate  thermal  collision  velocities 
of  4  kmy/s  (most  probable  collision  velocity  for  CO-O  at  Ttra  ~  13,000  K),  while  Figs.  1(b),  (d),  and  (f) 
show  the  RVEDF  at  a  hyperthermal  collision  velocity  of  10  km/s.  Since  the  current  work  focuses  on  assessing 
vibrational  relaxation  models,  we  sum  over  the  final  rotational  states.  For  reference,  the  energy  spacing 
between  the  first  excited  and  ground  vibrational  state  is  equivalent  to  the  translational  collision  energy 
between  CO  +  O  system  with  a  relative  velocity  of  about  2.2  km/s. 

Figure  2  shows  the  PVI-VEDF  at  at  a  relative  collision  velocity  of  4  km/s  that  was  calculated  by  summing 
over  all  possible  rotational  states  in  Figs.  1(a),  (c),  and  (e).  For  the  pure  excitation  case  (v{  =  0),  the 
available  total  energy  is  small  enough  to  restrict  the  VEDF  to  the  first  two  excited  states.  For  the  combined 
excitation/quenching  cases  (vi  =  4),  there  is  a  greater  probability  of  a  quenching  process  than  a  vibrationally 
excited  collision.  This  is  expected  because  the  initial  vibrational  energy  exceeds  60%  of  the  total  collisional 
energy.  This  VEDF  is  nearly  linear  in  vibrational  quanta.  Additionally  we  see  that  the  PVI-VEDF  remains 
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nearly  the  same  as  the  initial  rotational  temperature  is  changed.  This  gives  confidence  that  vibrational 
relaxation  models  may  not  require  other  internal  energy  states  as  inputs  to  maintain  physical  accuracy.  This 
is  the  case  for  both  vibrational  relaxation  models  studied  in  this  work. 

Figure  3  show  the  PVI-VEDF  at  a  relative  collision  velocity  of  10  km/s  that  was  calculated  by  integrating 
over  all  possible  rotational  states  in  Figs.  1(b),  (d),  and  (f).  At  this  collision  velocity,  we  see  the  pure 
excitation  case  (vi  =  0)  fills  many  of  the  excited  levels  with  measurable  cross-sections  up  to  Vf  =  17. 
Consistent  with  the  lower  relative  collision  velocity  case,  the  mixed  excitation/quenching  results  (vi  =  4) 
remain  similar  at  the  two  initial  rotational  temperatures.  However,  the  PVI-VEDF  is  split  nearly  equally 
between  excitation  and  quenching  with  the  jump  of  ±1  quanta  being  the  most  probable  final  state. 

III.  Vibrational  Relaxation  Models 

Modeling  of  collisional  vibrational  relaxation  within  the  DSMC  method  can  be  described  by  two  important 
steps  of  the  collision  process.  First,  a  decision  of  whether  vibrational  relaxation  should  occur  is  determined 
using  the  vibrational  cross  section,  which  may  be  a  function  of  many  collision  pair  parameters  such  initial 
vibrational  energy  content,  initial  rotation  energy  content,  and  relative  collision  energy.  This  vibrational 
cross  section  is  then  compared  to  the  total  cross  section,  atoti  to  form  a  probability  of  vibrational  energy 
exchange  that  is  used  within  a  Monte-Carlo  algorithm.  For  a  given  collision,  if  a  vibrational  energy  transfer 
process  is  selected,  then  post-exchange  vibrational  energy  states  are  selected  from  a  VEDF,  which,  in  this 
work,  is  termed  the  post  vibrationally-inelastic  VEDF  (PVI-VEDF),  which  again  can  be  dependent  on  a 
wide  range  of  parameters  such  as  initial  vibrational  energy  content  of  each  collision  pair,  initial  rotational 
energy  content  of  each  collision  pair,  and  relative  collision  energy.  Low-fidelity  models  typically  require  very 
little  information  to  predict  the  vibrationally  inelastic  cross  section  and  post  inelastic  vibrational  populations, 
while  high-fidelity  models  differentiate  many  collision  classes  through  state-resolved  vibrational  cross  sections 
and  PVI-VEDF.  In  this  description,  all  units  are  in  SI  unless  otherwise  noted. 

A.  Low-Fidelity  Vibrational  Relaxation  Models 

1.  Collision  Number  Selection 

Typically,  low-fidelity  vibrational  relaxation  models  form  the  probability  of  vibrationally  inelastic  cross 
section,  <ptotl  which  is  the  ratio  of  vibrational  cross  section  to  total  cross  section.  A  constant  value  based  on 
the  average  relaxation  time  can  be  applied  to  each  collision.  Alternatively,  the  probability  can  be  based  on 
a  mix  of  collisional  values  or  macroscopic  information  to  match  measured  relaxation  rates.  This  parameter 
can  be  freely  chosen  to  match  bulk  relaxation  and  has  little  effect  on  the  overall  model  form.  References  8 
and  9  provide  further  details  of  available  collision  number  models. 

2.  Post  Vibrationally  Inelastic  Population  Sampling 

After  a  vibrationally  inelastic  collision  is  selected,  energy  is  distributed  among  available  translational  and 
vibrational  energy  modes  using  a  modified  version  of  the  Larsen-Borgnakke2,3  (LB)  procedure  that  accounts 
for  quantized  vibrational  levels.  For  example,  the  LB  PVI-VEDF  for  a  harmonic  oscillator,  /  (£Vib)?  can  be 
defined  by  Eq.  1  where  eVIB  is  the  vibrational  energy,  Ec  is  the  total  collision  energy,  i  is  the  vibrational 
quantum  number,  k  is  the  Boltzmann  constant,  0V IB  is  the  characteristic  temperature  for  vibration,  and  S 
is  the  Dirac  delta  function.  This  model  includes  the  effect  of  the  initial  vibrational  energy  in  aggregate,  but 
cannot  differentiate  between  a  high  velocity  impact  collision  with  a  low  initial  vibrational  energy  from  a  low 
velocity  collision  with  a  high  initial  vibrational  energy  level. 

/  F  \Vtra- i 

/lb  (eVm)  ~  (  1  -  )  S  (eVIB  -  ikOy ib)  (1) 

B.  High-Fidelity  Vibrational  Relaxation  Model 

The  high-fidelity  vibrational  relaxation  model  uses  the  vibrationally  specific  transition  probabilities  of  the 
Forced  Harmonic  Oscillator  (FHO)  model  described  by  Adamovich  et  al.A  For  example,  the  probability  of 
vibrational  relaxation  for  an  atom-diatom  collision  is  described  by  Eq.  2  where  i  is  the  initial  vibrational 
state  of  the  diatom,  /  is  the  final  vibrational  state,  and  n  is  defined  as  minimum  of  i  and  /.  The  parameter 
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Figure  1:  Predicted  ro- vibrational  distribution  function  from  QCT  results  for  for  post  inelastic  collision 
states  given  the  relative  velocity  and  initial  vibrational  level  and  rotational  temperature  of  CO  molecule. 
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Figure  2:  Predicted  post-inelastic  vibrational  distribution  function  from  QCT  results  with  an  initial  relative 
velocity  of  4  km/s  and  specified  initial  vibrational  quantum  number  and  rotational  temperature  for  the  CO 
molecule.  Results  are  normalized  such  that  the  integral  of  the  VEDF  is  unity. 
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Figure  3:  Predicted  vibrational  distribution  function  from  QCT  results  with  an  initial  relative  velocity  of 
10  km/s  and  specified  initial  vibrational  quantum  number  and  rotational  temperature  for  the  CO  molecule. 
Results  are  normalized  such  that  the  integral  of  the  VEDF  is  unity. 
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£Col  is  defined  by  Eq.  3  where  SVT  is  a  steric  factor,  uj  is  the  oscillator  frequency,  rh  is  the  collision  reduced 
mass,  fi  is  the  oscillator  reduced  mass,  7  is  the  oscillator  mass  ratio,  h  is  the  Planck  constant,  v  is  the 
symmetrized  relative  velocity,  and  a  characterizes  the  intermolecular  potential. 


-Pvt  (i  /,  £col)  =  *!/!£co(  exp  (— eCOL) 


(-l)r 


fr'o  P  (*  -  rV-  (/  -  r)! 


(2) 


_  c  47r3uj  (™2/v)  7  firu\  ^ 

^col  —  *^vt  9,  smh  (  _  j  (3) 

azri  \avJ 

For  diatom-diatom  collisions,  the  probability  of  vibrational  relaxation  is  described  by  Eq.  4  using  the 
same  definitions  as  the  atom-diatom  collision  and  Cca  b  is  the  transformation  matrix10  whose  coefficients  are 
defined  by  Eq.  5.  The  parameter  pCOL  is  defined  by  Eq.  6  where  Svv  is  a  steric  factor,  and  ujn  is  the 
oscillator  frequency  of  molecule  n. 


Pwt  [}!•)  ^2  ~ ^  /l?  f2i  ^COL)  Pool)  — 


E”=0  (-  !)ii+i2  _r  cfj+fi 


yr+l,22  +  lwr+l,/2+l 


exp  (  i/>coL  (/l  f 2  r))  Pvt  (H  4“  ^2  r  *  fl  P  /2  ^5  2£COL) 


(4) 


p  1 

V2 

r  1 

V2  j- 1 

p  1 

1- 

Ckj  =  2“n/2 

n 

k-  1 

n 

_  J  ~  1  _ 

E(-1)“^r 

r=0 

n  —  k  +  1 
.  4  -  r  ~  1 

k  —  1 

r 

(5) 


/  a2i;2  \ /2 

PcOL=  (Sw -  (6) 

V  ^1^2  ) 

In  the  FHO  model,  the  vibrational  cross-sections  and  the  PVI-VEDF  are  inherently  linked.  The  steric 
factors,  SVT  and  Svv ,  and  the  parameter  a  determine  the  rate  of  vibrational  relaxation,  while  the  inherent 
form  of  the  transition  probability  determines  the  PVI-VEDF.  Because  of  the  complexity  of  the  mathematical 
form,  Boyd  and  Josyula1  suggested  accumulating  the  transition  probabilities  into  two  tables.  For  example, 
for  an  atom-diatom  collision  given  an  initial  state,  (i,  £Col)>  the  probability  of  any  change  in  vibrational 
quantum,  PVt-tol,  can  be  accumulated  into  the  first  table  by  summing  over  all  possible  final  states  as  seen 
in  Eq.  7. 


f max 

Pw T  — TOL  =  ^  ^  P'V'T  ^  />£ COl) 


(7) 


Then,  a  second  table  gives  the  PVI-VEDF,  /PVi-vedf  (i  /,£Col),  for  each  initial  state  as  seen  in  Eq. 

-Pvt  ^  f  5  ^COL)  i  7^  /) 


.  ft  ^col) 


(8) 


This  lends  itself  to  a  similar  computational  algorithm  used  in  the  low-fidelity  vibrational  relaxation  models 
where  the  first  table  provides  a  vibrational  collision  number  that  is  dependent  on  initial  conditions.  If  selected 
for  a  vibrationally  inelastic  collision,  the  second  table  provides  the  scaled  PVI-VEDF  that  is  sampled. 


IV.  Assessment  of  Models 

In  this  section,  we  compare  the  vibrational  relaxation  models  to  the  quasi-classical  trajectory  (QCT) 
calculations  of  Braunstein  and  Duff6  for  the  C0(^,TROt  =  300 K)  +  O  — >>  CO (vf,jf)  +  O  system  over  a 
wide  range  of  relative  collision  velocities.  Both  models  have  internal  parameters  that  can  be  set  to  match 
consistent  vibrational  collision  numbers.  However,  the  vibrational  relaxation  model  form  dictates  the  PVI- 
VEDF  and  therefore  the  evolution  of  the  detailed  vibrational  relaxation  process.  Then,  we  isolate  the  effects 
of  the  model  form  for  a  zero-dimensional  vibrational  relaxation  case  and  compare  the  evolution  of  the  VEDF 
and  radiative  emission  predicted  by  each  model.  Finally,  we  present  how  these  differences  affect  the  predicted 
radiative  emission  spectrum. 
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A.  Initial  Model  Comparison 

Figures  4  and  5  compare  the  PVI-VEDF  predicted  by  the  FHO  and  LB  with  the  QCT  data  at  an  initial 
rotational  temperature  of  300  K;  initial  vibrational  levels  of  0,  2,  and  4;  and  relative  collision  velocities  of 
4,  6,  8,  and  10  km/s.  In  general,  the  LB  model  predicts  that  the  probability  for  the  final  vibrational  level  is 
inversely  proportional  to  the  final  vibrational  level.  In  contrast,  the  FHO  model  generally  predicts  that  the 
probability  for  each  final  state  is  inversely  proportional  with  the  vibrational  quanta  jump.  In  general,  the 
LB  model  is  in  better  agreement  with  the  QCT  data  at  lower  collision  energies,  while  the  FHO  model  is  in 
better  agreement  with  the  QCT  data  at  higher  collision  velocities. 

Figure  4(a)  shows  the  PVI-VEDF  from  an  initial  ground  vibrational  state  and  collision  velocity  of  4  km/s. 
Here,  both  models  and  QCT  results  predict  that  almost  all  molecules  are  excited  to  the  first  two  excited 
states,  with  a  single  quantum  excitation  having  a  probability  greater  than  85%.  However,  compared  to  the 
QCT  results,  both  vibrational  relaxation  models  under-predict  the  probability  of  a  double  quanta  excitation 
and  consequently  over-predict  a  single  quanta  excitation.  Figures  4(b)  and  5(a)  and  (b)  show  similar  trends 
with  an  increase  in  the  probability  of  multiple  quantum  jumps  that  increases  with  the  increase  in  relative 
collision  velocity.  The  FHO  and  LB  tend  to  predict  fewer  multiple  state  transition  than  the  QCT  data,  with 
the  LB  model  maintaining  a  better  level  of  agreement  with  the  QCT  data. 

Figure  4(c)  shows  the  PVI-VEDF  for  a  CO  molecule  initially  in  the  second  excited  vibrational  state  and 
collision  velocity  of  4  km/s.  Both  the  LB  model  and  QCT  data  predict  a  probability  that  is  proportional  to 
the  final  state,  while  the  FHO  model  predicts  a  probability  that  is  proportional  to  the  quanta  jump.  Figures 
4(d)  and  5(c)  and  (d)  show  similar  trends  for  both  LB  and  FHO  models  with  an  increase  in  possible  final 
vibrational  states  with  an  increase  in  relative  collision  velocity.  Interestingly  at  the  higher  velocities  of  8  and 
10  km/s,  the  QCT  data  shows  a  change  in  trend  where  the  probability  of  a  post-collision  state  in  the  ground 
vibrational  state  is  less  than  the  first  and  third  excited  state.  This  results  in  the  FHO  model  qualitatively 
showing  the  same  trend  as  the  QCT  data.  With  the  increase  in  available  post-collision  states,  the  FHO 
model  begins  to  show  local  maxima  that  are  away  from  the  single  quanta  jump.  Similar  local  maxima  have 
been  calculated  with  the  FHO  model  near  resonant  collision  velocities.11 

Figures  4(e)  and  (f)  show  the  PVI-VEDF  from  Vi  =  4  and  collision  velocity  of  4  and  6  km/s,  respectively. 
The  QCT  data,  LB,  and  FHO  models  all  predict  similar  trends  as  they  do  for  the  collision  process  from 
Vi  =  2.  Figures  5(e)  and  (f)  show  the  PVI-VEDF  from  the  initial  fourth  excited  vibrational  state  with 
collision  velocities  of  8  and  10  km/s,  respectively.  Again,  we  see  very  similar  trends  from  the  second  excited 
vibrational  state.  Now,  the  FHO  model  predicts  more  pronounced  local  maxima  away  from  the  single  quanta 
jump.  We  also  see  that  these  maxima  can  shift  to  higher  states  as  the  relative  collision  velocity  increases.  For 
quenching  processes,  the  QCT  data  predicts  slight  maxima,  however  the  QCT  data  always  monotonically 
decreases  as  the  vibrational  energy  level  jump  increases  for  all  excitation  processes.  As  with  all  other  collision 
classes  studied,  the  LB  model  consistently  predicts  that  the  probability  is  proportional  to  the  final  vibrational 
state. 

B.  Vibrational  Relaxation  Comparison 

The  effects  of  the  model  form  can  be  isolated  from  the  collision  number  using  numerical  simulations.  Both  the 
LB  and  FHO  model  are  used  within  the  DSMC  method  to  predict  adiabatic  relaxation.  A  zero  dimensional, 
unsteady  configuration  in  which  the  energy  modes  of  the  gas  are  initialized  to  Boltzmann  distributions 
specified  by  a  mode  temperature  is  used  to  model  the  basic  vibrational  relaxation  process.  The  initial 
conditions  for  the  adiabatic  relaxation  case  presented  here  include  an  initial  pressure  of  p  =  101,325  Pa, 
translational  temperature  of  Ttra  =  5,  000  K  and  vibrational  temperature  of  TVIB  =  1, 000  K.  The  gas  species 
concentration  consists  of  50%  CO  and  50%  O  by  mole  fraction.  Vibrational  relaxation  between  CO  and  CO 
is  suppressed  to  isolate  the  effects  of  the  CO  +  O  collisions.  In  addition,  the  effects  of  radiative  decay  on 
the  VEDF  is  ignored.  Both  simulations  include  consistent  rotational  relaxation  models  for  the  CO  +  O  and 
CO  +  CO  collision  classes. 

Figures  6  and  8  show  the  time  evolution  of  the  translational  and  vibrational  temperature  obtained  with 
DSMC  using  both  the  FHO  and  LB  models.  Here,  the  translational  temperature  is  calculated  using  thermal 
velocity  components  as  seen  in  Eq.  9  where  c  is  the  particle  velocity,  m  is  the  particle  mass,  and  k  is  the 
Boltzmann  constant.  The  vibrational  temperature  is  calculated  using  the  average  vibrational  energy  per 
molecule,  e^,  as  shown  in  Eq.  10  where  gi  is  the  degeneracy,  and  is  the  vibrational  energy  of  level 

i.  For  this  study,  the  collision  number  for  the  LB  model  are  selected  to  match  the  macroscopic  relaxation 
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^  final  ^  final 

(e)  Vi  =  4,  grel  =  4  km/s  (f)  V|  *=  4,  grel  =  6  km/s 

Figure  4:  Predicted  PVI-VEDF  from  vibrational  relaxation  models  and  QCT  results  with  an  initial  rotational 
temperature  of  300  K;  initial  vibrational  level  of  0  (top),  2  (middle)  and  4  (bottom);  and  relative  velocity  of 
4  km/s  (left)  and  6  km/s  (right). 
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final 


(a)  vi  =  0,  grei  =  8  km/s 


^final 

(b)  vi  =  0,  grei  =  10  km/s 


0.30 


^final 

(c)  Vi  =  2,  grei  =  8  km/s 


^  final  ^final 

(e)  Vi  =  4,  grei  =  8  km/s  (f)  v.  =  4,  grel  =  10  km/s 

Figure  5:  Predicted  PVI-VEDF  from  vibrational  relaxation  models  and  QCT  results  with  an  initial  rotational 
temperature  of  300  K;  initial  vibrational  level  of  0  (top),  2  (middle)  and  4  (bottom);  and  relative  velocity  of 
8  km/s  (left)  and  10  km/s  (right). 


9  of  14 


American  Institute  of  Aeronautics  and  Astronautics 
Distribution  A:  Approved  for  public  release;  distribution  unlimited. 


predicted  by  the  FHO  model.  Indeed,  Figure  6  shows  that  careful  collision  number  selection  can  match  the 
macroscopic  temperatures  predicted  using  each  model  over  the  entire  relaxation  process.  It  is  important 
to  note,  that  although  the  time-scales  in  the  selected  problem  are  too  short  for  available  measurement 
techniques,  these  time-scales  increase  through  a  decrease  of  the  initial  pressure,  and  therefore,  increase  in 
the  collision  time-scale. 


Ttra  -  3k  [c2 


C21 


&vib  — 


Tn=ieU  (H^vibA  exP  Tyjb) 

EZTls  9i  exP  {-e'’ib’i/kTvib) 


(9) 

(10) 


Figure  6:  Evolution  of  adiabatic  vibrational  relaxation  (excitation)  using  LB  and  FHO  vibrational  models 
within  the  DSMC  method. 

Despite  the  excellent  agreement  in  the  macroscopic  temperatures,  the  form  dictated  by  the  FHO  and 
LB  models  result  in  a  pronounced  difference  in  the  prediction  of  the  vibrational  energy  distribution  function 
(VEDF).  For  example,  Figures  7  and  8  (left)  show  the  evolution  of  the  VEDF  using  each  model.  For 
reference,  the  equilibrium  Boltzmann  VEDF  is  shown  at  the  bulk  vibrational  temperature.  As  the  VEDF 
transitions  from  the  low  to  high  temperature,  both  the  FHO  and  LB  models  predict  a  bi-modal  distribution 
with  a  large  percentage  of  low  energy  particles  and  a  high  energy  tail.  Compared  to  the  FHO  model,  the 
LB  model  predicts  a  larger  change  in  slope  of  the  VEDF  that  can  described  as  the  mode  temperature.  This 
is  due  to  a  larger  disparity  in  the  mode  temperatures  using  the  LB  model  and  results  in  the  FHO  model 
predicting  higher  populations  of  the  ground,  first,  and  second  excited  states,  while  the  LB  model  predicts 
higher  populations  in  the  third  and  higher  excited  states.  Both  models  reach  an  equilibrium  VEDF  that 
can  be  described  by  a  single  temperature.  This  apparent  bi-modal  VEDF  can  be  characterized  by  fitting 
each  VEDF  to  the  functional  form  shown  in  Eq.  11  where  two  temperatures,  i  and  2,  and  relative 
populations,  niand  77,2,  are  selected  to  fit  each  sampled  distribution  function,  /Vib-sample-  These  fits  are 
also  shown  in  Figures  7  and  8(left). 

~Evib  ~Evib 

fv IB-SAMPLE  ~  /v IB-FIT  ^  nl  6  kTvib’ 1  +  Tl2^'hT’V%h'1  (H) 

Tables  1  and  2  list  the  coefficients  for  the  VEDF  fits  at  each  time-sample  slice  for  the  LB  and  FHO  models, 
respectively.  The  LB  model  predicts  a  VEDF  that  consists  of  two  populations  that  are  at  the  initial  and  final 
temperatures  with  a  steady  depletion  of  one  to  the  other.  Even  at  short  times,  the  high  temperature  tail  can 
be  described  with  a  temperature  that  is  within  10%  of  the  final  temperature.  In  contrast,  the  FHO  model 
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predicts  a  VEDF  that  is  more  gradual  in  the  transition  of  modes  with  both  mode  temperatures  increasing 
as  the  bulk  vibrational  temperature  approaches  equilibrium. 

The  distribution  of  vibrational  energy  is  used  to  predict  the  spectral  emission.  The  fit  parameters  are 
used  to  construct  a  gas  mixture  for  a  line  of  sight  that  runs  through  a  uniform  gas.  These  lines  of  sight  are 
provided  to  an  approximate  radiative  transfer  algorithm  to  predict  IR  emission  for  fundamental  and  hot  band 
emission  at  moderate  spectral  resolution.12  This  work  assumes  that  the  gas  is  optically  thin.  The  evolution 
of  the  spectral  emission  for  the  fundamental  and  first  overtone  can  be  seen  in  Figures  7  and  8 (right).  Even  at 
the  initial  sample  slice  of  1  x  10-10  s,  the  predicted  emissive  spectrum  from  FHO  and  LB  models  differ.  The 
differences  can  be  traced  back  to  the  VEDF  which  is  shown  to  the  left  of  each  spectrum.  For  example,  the 
FHO  VEDF  results  in  a  higher  spectrum  across  most  of  the  fundamental  transition  states  1,500-2,300 
cm-1)  where  the  1  — 0  emission  dominates  the  spectra.  In  this  region,  the  peaks  for  all  upper  vibrational 
transitions,  e.g.  2  — >•  1;  3  — >•  2;  etc.,  are  combined  with  the  P-branch  of  the  1  — >•  0  transition.  The  first 
overtone  3, 000  —  4, 400  cm-1)  clearly  shows  the  bandhead  structure.  This  structure  shows  clear  emission 
peaks  for  the  2  — >>  0  4,  370  cm-1),  3  — 1  (~  4, 315  cm-1),  4  — >>  2  4,  260  cm-1),  etc  transitions.  As  the 

LB  VEDF  predicts  higher  populations  for  most  excited  vibrational  states,  it  predicts  higher  radiance  than 
the  FHO  model  for  most  of  the  overtone  emission,  outside  of  the  2—^0  peak.  This  is  expected  since  the 
FHO  VEDF  predicts  a  higher  population  for  the  second  excited  vibrational  state. 

These  trends  continue  throughout  the  relaxation  process.  Until  both  VEDFs  reach  the  Boltzmann  distri¬ 
bution  function  at  equilibrium,  the  FHO  always  predicts  a  higher  population  in  the  first  excited  vibrational 
state,  which  results  in  the  FHO  predicting  a  higher  radiance  over  most  of  the  fundamental  band.  The  LB 
VEDF  predicts  higher  populations  for  all  most  excited  states  which  results  in  a  higher  emission  over  most 
of  the  overtone  band.  Although  both  models  differ  from  the  Boltzmann  VEDF  at  the  bulk  vibrational  tem¬ 
perature,  the  LB  model  predicts  a  higher  level  of  bi-modal  behavior  throughout  the  relaxation  process.  At 
equilibrium,  both  models  predict  the  same  Boltzmann  VEDF  which  results  in  a  very  high  level  of  agreement 
in  the  predict  spectra. 


Table  1:  Bi-Modal  VEDF  LB  Fit  Coefficients 


Time  [s] 

711  f  riTOT 

TVib.I  [K] 

n<1  /  riTOT 

TVIB,2  [K] 

l  x  icr10 

.96 

1009 

.04 

3981 

1  x  10"9 

.88 

1029 

.12 

4001 

3  x  10"9 

.65 

1117 

.35 

4100 

i  x  icr8 

.22 

1249 

.78 

4212 

3  x  icr8 

1 

4288 

3  x  lO"7 

1 

4300 

Table  2:  Bi-Modal  VEDF  FHO  Fit  Coefficients 


Time  [s] 

ni/ riT  OT 

Tvib.I  [K] 

n2/ TlTOT 

Tvib,2  [K] 

1  x  10-10 

.77 

1018 

.23 

2064 

1  x  10"9 

.60 

1201 

.40 

2476 

3  x  10“9 

.19 

1418 

.81 

3089 

1  x  10“8 

.05 

1818 

.95 

4010 

3  x  10“8 

.01 

155 

.99 

4362 

3  x  10-7 

.02 

145 

.98 

4396 

V.  Summary  and  Conclusions 

The  form  of  two  vibrational  relaxation  models  that  are  commonly  used  in  DSMC  and  CFD  simula¬ 
tions  have  been  compared.  In  particular,  this  work  focused  on  examining  the  effect  of  the  model  form  for 
both  the  Larsen-Borgnakke  (LB)  and  Force  Harmonic  Oscillator  (FHO)  vibrational  relaxation  models  and 
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Probability  Density  Probability  Density  Probability  Density 


(a)  VEDF,  t  =  1  x  10  10  s  (b)  Spectra,  t  =  1  x  10  10  s 


(c)  VEDF,  t  =  1  x  10  9  s  (d)  Spectra,  t  =  1  x  10  9  s 


(e)  VEDF,  t  =  3  x  10  9  s  (f)  Spectra,  t  =  3  x  10  9  s 

Figure  7:  Evolution  of  the  VEDF  (left)  and  radiative  emission  (right)  using  FHO  and  LB  vibrational  models 
for  adiabatic  relaxation. 
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(a)  VEDF,  t  =  1  x  10  8  s  (b)  Spectra,  t  =  1  x  10  8  s 
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(d)  Spectra,  t  =  3  x  10  8  s 


5000 


(e)  VEDF,  t  =  3  x  10  7  s  (f)  Spectra,  t  =  3  x  10  7  s 

Figure  8:  Evolution  of  the  VEDF  (left)  and  radiative  emission  (right)  using  FHO  and  LB  vibrational  models 
for  adiabatic  relaxation  (continued). 
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their  effect  on  the  vibrational  energy  distribution  function  (VEDF)  and  radiative  emission  spectra.  Avail¬ 
able  ro-vibrationally  resolved  cross  sections  from  state  of  the  art  QCT  calculations  were  presented  for  the 
CO Trot  =  300  K)  +  O  -A  CO (vf,jf)  +  O  system.  This  data  was  reduced  to  post-vibrationally  inelastic 
VEDF  (P VI- VEDF)  for  direct  comparison  with  each  model  predictions  over  a  wide  range  of  collision  energies. 
Comparison  between  each  model  with  the  QCT  data  showed  that  the  LB  model  is  in  better  agreement  with 
the  QCT  data  for  lower  collision  velocities,  while  the  FHO  model  is  in  better  agreement  for  high  collision 
velocities.  Each  model  was  used  within  the  DSMC  method  for  an  adiabatic  thermal  relaxation  simulation. 
Even  when  the  collision  numbers  were  selected  to  predict  the  same  bulk  temperature  relaxation  profile,  the 
LB  and  FHO  models  predicted  very  different  VEDFs.  Both  models  predicted  bi-modal  VEDFs,  however  the 
evolution  of  the  bi-modal  temperatures  and  populations  were  different  due  to  the  model  form.  The  disparity 
in  the  predicted  VEDF  resulted  in  differences  in  the  predicted  radiative  emission  spectra. 

Future  work  will  aim  to  extend  the  evaluation  of  each  model  with  additional  available  high  fidelity 
chemistry  data  and  study  the  effect  of  each  model  on  a  wide  range  of  compression  and  expansion  flows.  By 
leveraging  recent  developments  that  use  classical  trajectory  calculations  directly  within  the  DSMC  collision 
kernel,13  the  vibrational  relaxation  models  can  be  compared  to  QCT  predictions  for  problems  similar  to  the 
test  case  used  in  this  paper.  In  addition,  study  of  the  effect  of  the  model  form  on  chemistry  would  provide 
further  understanding  of  coupled  vibration-dissociation  processes  that  are  important  in  many  high  altitude 
flows.  Finally,  further  development  of  vibrational  relaxation  models  depends  on  high  quality  state-resolved 
experimental  measured  data.  These  initial  results  have  shown  that  each  model  results  in  significant  differences 
in  a  measurable  optical  quantity.  We  note  that,  although  the  time-scales  for  the  relaxation  process  in  this 
study  are  quite  small,  they  can  be  increased  through  a  reduction  in  pressure  and  the  subsequent  increasein 
the  collision  time-scale.  Spectrally-resolved  measurements  for  similar  cases  under  non-equilibrium  conditions, 
such  as  in  shock  tube  experiments  can  provide  further  insight  into  inelastic  collision  processes. 
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